Singular Dynamics of Underscreened Magnetic Impurity Models 



O 

o 

> 

O 

O 

in 



Winfried KollerQ Alex C. HewsonQ and Dietrich MeyeJl 

Department of Mathematics, Imperial College, London SW7 2BZ, United Kingdom 

We give a comprehensive analysis of the singular dynamics and of the low-energy fixed point of 
one-channel impurity s-d models with ferromagnetic and underscreened antiferromagnetic couplings. 
We use the numerical renormalization group (NRG) to perform calculations at T — 0. The spectral 
densities of the one-electron Green's functions and t-matrices are found to have very sharp cusps at 
the Fermi level (w = 0), but do not diverge. The approach of the Fermi level is governed by terms 



proportional to 1/ln (uj/Tq) as 



0. The scaled NRG energy levels show a slow convergence 



as 1/(JV + C) to their fixed point values, where N is the iteration number and C is a constant 
dependent on the coupling J from which the low energy scale To can be deduced. We calculate also 
the dynamical spin susceptibility, and the elastic and inelastic scattering cross-sections as a function 
of lo. The inelastic scattering goes to zero as u> — > 0, as expected for a Fermi liquid, but anomalously 
slowly compared to the fully screened case. We obtain the asymptotic forms for the phase shifts for 
elastic scattering of the quasiparticles in the high-spin and low-spin channels. 
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PACS numbers: 75.20.Hr,71.27.+a 



I. INTRODUCTION 

The breakdown of Fermi liquid behaviour in the neigh- 
bourhood of a quantum critical point (QCP) has been 
observed in an increasing number of heavy fermion ma- 
terials in recent yearsi, but so far has not received an 
adequate explanation. The QCP is the point at which 
the transition temperature Tn for long range magnetic 
order, usually antiferromagnetic, is such Tn = 0, having 
been driven down via the application of pressure or alloy- 
ing to a magnetically ordered material (see for example 
Refs. yiij). In these circumstances the Wilson theory for 
critical behaviour has to be extended to include quan- 
tum as well as thermal fluctuations of the order param- 
eter. Such generalizations have been carried out^^ but 
the results do not appear to be consistent with the ex- 
perimental observations. The experimental evidence, for 
instance from neutron scattering in the heavy fermion 
alloy CeCu5.9Auo.1- where ui/T scaling is found at the 
QCP, indicates that the critical behaviour is associated 
with very short range fluctuations in contrast with crit- 
ical phenomena driven by purely thermal fluctuations, 
where the shorter range fluctuations are only important 
in so far as they modify the longer range fluctuations. 

One conjecture that has been put forward as a basis 
for a possible theoretical explanation is that the local 
moments, free from the constraints of the magnetic order, 
yet not fully screened as in a conventional S — 1/2 Kondo 
model, cause singular scattering leading to a breakdown 
of Fermi liquid theoryi Such a breakdown can occur 
in impurity models with a pseudogap at the Fermi level. 
This type of model has been extensively studied as a 
model which has a local QCP, see e.g. Refs. IllTnllTiI 
For such a model to explain the observed critical 
behaviour, however, an explanation would be required 
for the origin of the pseudogap. 

The pseudogap model is not the only form of impu- 
rity model, however, that can lead to singular scatter- 
ing. It was recognized by Coleman and Pepm—i that 



the underscreened s-d model has free spins at T = 
that can give rise to singular scattering. The s-d 
model is underscreened when 2S, where S is the quan- 
tum number of the antiferromagnetically coupled local 
spin, exceeds the number of conduction electron chan- 
nels. The thermodynamical properties of this model have 
been known for over 20 years from exact Bethe ansatz 
calculationsi^i^iSiii. The dynamical behaviour, how- 
ever, has only been studied recently and mainly via gen- 
eralizations to SU(N) versions of the model in the large 
N limit, where mean field methods can be applied to- 
gether with estimates of the leading Gaussian correction 
terms 13,18,19 . There has also been some recent NRG work 
and calculations of the spinon density of states via the 
Bethe ansata2&. 

The underscreened s-d model might be realizable for 
certain quantum dot configurations, as discussed in 
Ref. This possibility has already been considered in 
a recent preprint, Ref. l22t 

Previous NRG calculations for the underscreened 
models 20, — have been confined to the calculation of the 
energy levels, using the original Wilson procedural. 
Here we apply the generalization of NRG approach^^ to 
carry out extensive calculations of the dynamical quanti- 
ties for the one-channel underscreened model for various 
values of the spin S. We also discuss the model with a 
ferromagnetically coupled spin for comparison. The cal- 
culations are for the Hamiltonian of the s-d model in the 
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interacting with a localized spin S, 

H sd = ^2 Jkk' (S + c k ^c k ,^ + S c\^c k , ^ 



kk' 



(1) 



(2) 



We will assume a separable form for the interaction, 
Jkk' = Joikotk' , and take d a = J2k a k c k a t° be the site of 
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FIG. 1: This figure illustrates the notation used in this paper 
for the s-d- model (top) and the Anderson model (bottom). 
The site coupling to the impurity is denoted by d, and in the 
case of the Anderson model the Coulomb repulsion U acts at 
the site / only. 

the lattice that couples to the impurity spin, i.e., the first 
site of a conduction electron chain as used in an NRG 
calculation. For S — 1/2 and J > this is commonly 
known as the Kondo Hamiltonian. We use a flat density 
of states po = 1/(2D) in the interval [-D, D] with D = 1 
and concentrate on the particle-hole symmetric case only. 
The notation used in this paper is illustrated in Fig. ^ 
It should be noted that the impurity spin is located at 
site N = -2. 

We briefly outline the structure of the paper: In section 
|n|we analyse the renormalization group flows of the low- 
lying excitations in terms of an effective exchange model, 
and determine the renormalized coupling in terms of an 
effective energy scale To- We then calculate the spectral 
densities of the local one-particle Green's functions and 
t-matrices for T = as a function of u> in section ITTT1 the 
elastic and inelastic scattering cross-sections in section 
IIVI and the dynamical spin susceptibilities in section \V\ 
In the final section we review the overall picture that 
these results give into the physics of this class of impurity 
models, and discuss the relation between our conclusions 
and previous work on this topic. 



II. FIXED POINT ANALYSIS 

As a preliminary to the calculation of the dynamics, 
we look at the NRG flows and the approach to the fixed 
point. This will clarify the nature of the fixed point and 
the form of the leading correction terms. For this analysis 
we use a relatively large value of the NRG discretisation 
parameter A = 3.5 and retain 600 states. Due to the slow 
convergence to the fixed point, we find it necessary to do 
up to 250 iterations. 



A. Ferromagnetic coupling (J < 0) 

We begin the fixed-point analysis with the ferromag- 
netic case. Here 2 ^ the fixed-point Hamiltonian corre- 
sponds to free fermions on the uncoupled chain (J = 0) 
in the limit N — > oo with the additional degeneracy due 
to the uncoupled spin S. The Hamiltonian for the uncou- 
pled chain for finite N can be diagonalized and written 



in the form 

(N+2)/2 

A -(N-i)/2 J- (E^kN)pt, a P r ^+E^(N)ht t(r h rta ), 

r=l 

(3) 

where pj, a , p ra , and h' Tcr , h ra , are the creation and an- 
nihilation operators for the single particle and hole ex- 
citations, and A-W-WEjftiN) and A~( N ~ 1 V 2 Ej®.(N) 
are the corresponding excitation energies relative to the 
ground or vacuum state |0); the scale factor A - ^ - *" 2 
is due to the fact that the energies are calculated for 
the rescaled Hamiltonian. We have assumed particle- 
hole symmetry, Ep°}(N) — E^l(N), so we can drop the 
p and h labels in denoting the levels. The levels are or- 
dered such that £y+i(A0 > e1 0) (N) and the scaling is 
such that for r <C N E^ (N) is of order 1 and has a finite 
value E* in the limit N — > oo. 

The leading correction terms to the fixed point are 
expected to be the simplest local terms consistent with 
the symmetry of the model (see Wilson's original pa- 
per— for a more detailed discussion). In this case they 
are expected to be of the same form as the original 
Hamiltonian 23 , and for the N + 2 site system can be 
expressed in the form, 

H C (N) = J(N)(S + d\d 1 + S-d\d l + S z (d]d^ - djd^)). 

(4) 

An on-site potential scattering term is excluded due to 
particle- hole symmetry. As H C (N) describes the excita- 
tions relative to the ground state of the uncoupled chain 
(iV + 2 sites) , it has to be normal ordered and written in 
terms of the single particle and hole operators, using 

4 = a 0r(p r t a + h r<a ) (5) 

r 

where the coefficients ao r will depend upon N. 

For the particle excitations only, H C (N) takes the form 

H C {N) =J{N) a 0r a 0rl (S+pl lPr ,^ + S-pl j1 p r , d 

r,r' 

+ S z {pl A p r , A -pt tl Pr',l)) ■ (6) 

There is an additional similar term for the hole excita- 
tions, and another one containing combinations of par- 
ticle and hole excitations. We will concentrate on the 
effects of these terms on the lowest single-particle and 
two-particle excitations of the N + 2 site system. 

In the ferromagnetic case, we consider even iterations 
in which case the ground state has total charge Q — 
and total spin 5tot = S. The lowest single particle 
excitations can have quantum numbers 5tot = S ±1/2 
and Q = 1. The corresponding energies are denoted by 
E Pi i(N, S ± 1/2). As these terms become asymptotically 
small in the limit N — > oo we calculate the shift in this 
excitation to first order in J(N). The total energy of the 
scaled excitation with spin S + 1/2 is then given by 

E pA (N,S+ 1/2) = (AO + /(AOSc&A^- 1 )/ 2 (7) 
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and that for spin 5—1/2 

E pA (N,S- 1/2) = E^l(N) - J(N)(S + l^A^ 1 )/ 2 . 

(8) 

From the fact that the high-spin excitation is found to lie 
lowest we infer that the effective interaction J(N) must 
be ferromagnetic (J(N) < 0). The energy difference be- 
tween high-spin and low-spin excitations is given by 



E pA (N,S + 1/2) 



E P ,i(N,S- 
: J(N)(2S- 



-1/2) 

lXA^- 1 )/ 2 



(9) 



This difference approaches its fixed point value, which 
is zero, very slowly, whereas the factor a§ 1 A^ Ar_1 - )//2 
rapidly reaches a finite fixed point value denoted by 
a 2 ,! = limAr^ooa^A^- 1 )/ 2 « 0.27 for A = 3.5. This 
suggests that the asymptotic form of J(N) for large N is 
given by 



J(N) = 



A 



N + C(J) 



(10) 



where A > is a constant. The slow fall off of J(N) with 
N, J(N) oc 1/N, indicates that the exchange term is a 
marginally irrelevant interaction 2 ^. 

We can translate the ./V-dependence of the effective 
parameter J{N) into a frequency or temperature scale 
using cjjv = r/DA^^^ 1 ^ 2 or T/v = ujn, where rj is an 
appropriately chosen constant of order unit y 24 ' 2 ^ which 
we take as n = 1. The relation is 



N 



In A 



In 



nD 



(11) 



Upon inserting this into Eq. (fTUf) and expressing C(J) in 
terms of an energy scale Tq, we obtain 
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(12) 



- ln(w/»jD) + HTo/vD) ' 

with A = A In A/2 and 

T a = nDe- c ^/ 2lnA . (13) 

Thus we find J{uo) oc l/m(w/T ) or J(T) cx l/ln(T/T ). 
In Sec. E| we show that this result is in line with the 
fitting to the low frequency part of the spectrum of the 
dynamic susceptibility for T = 0. It is also in line with 
the exact low temperature thermodynamics known from 
the Bethe ansata 14 i 15 4 6 . 

For small values of J, the dominant contribution to 
C(J) is proportional to 1/J. Hence we rewrite it in the 
form 



C(J) = 



B(J) 
J 



(14) 



with the parameter B(J) only weakly dependent on J. 
The energy scale T reads 



T Q = r]De 



B(J)/(2J) InA 



(15) 
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FIG. 2: The flow of the energy difference TV x (E P>1 (S+1 /2) - 
E Pt i(S — 1/2)) of high-spin and low-spin one-particle exci- 
tations for different values of S and ferromagnetic coupling 
J = -0.1. The symbols (x for S = 1, o for S = 3/2, and □ 
for 5 = 2) show the flows rescaled by l/(S + 1/2). 
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S=l/2 


S=l 


S=3/2 


S=2 


-0.005 


1.5990 


1.5991 


1.5992 


1.5994 


-0.010 


1.6034 


1.6036 


1.6039 


1.6043 


-0.025 


1.6226 


1.6223 


1.6219 


1.6214 


-0.050 


1.6637 


1.6586 


1.6516 


1.6426 


-0.100 


1.7589 


1.7270 


1.6827 


1.6261 


-0.200 


1.9540 


1.7970 


1.5789 


1.3017 



TABLE I: Values for B(J) obtained by numerical fitting of 
the low-energy behaviour (iterations 50 — 250) to Eq. (1101) and 
Eq. I|14|l using 250 iterations at A = 3.5. 



Figure |21 shows the flows of the energy difference N x 
{E P {S + 1/2) - E p (S - 1/2)) for J = -0.1 and different 
values of S. The rescaled curves (symbols) show that A 
and B(J) are indeed independent of S. We can obtain 
the numerical values for A and B{ J) by inserting Eq. I|14|) 
into Eq. lfTU|) and the result into Eq. ©. Fitting this to 
the energy difference obtained from the NRG, we obtain 
a value of A w 1.56 ± 0.01 independent of J and S for a 
small J. The values for B(J) are given in Table ||| 

Our value of A sa 1.56 yields an A = 0.977 w 1. 
Eq. iTHH) thus simplifies to 



J(w) = l/ln(w/T ), 



(16) 



which agrees with the scaling result used for the spin 
susceptibility in Eq. I|C5(I . 

In the limit N — * oo, the energy E pPy i(N, S) of the 
lowest two-particle excitation with quantum numbers 
Q = 2 and S to t = S is equal to 2 x E p °l(N) of the non- 
interacting case. The difference E pp ^(N, S) — 2E^\(N) 
falls off relatively slowly with N, approximately as 1/N 2 , 
which we can interpret as a second order effect of J(N). 



50 



100 



150 



200 



250 



-3.5 



50 



100 



150 



200 



250 



FIG. 3: The flow of the energy difference N x (E pA (S) - 
E P: i(S — 1)) of high-spin and low-spin one-particle excitations 
for different values of S and antiferromagnetic coupling J = 
+0.2. The symbols (x for S = 3/2, o for S = 2, and □ for 
S — 5/2) show the flows rescaled by 1/S. 



FIG. 4: The flow of the energy difference N x (E p> i(S) - 
E Pt i(S— 1)) of high-spin and low-spin one-particle excitations 
for 5 = 3/2 different values of the antiferromagnetic coupling 
J. 



Following the analysis used by Hofstetter and Zarand 27 
we can deduce the asymptotic form the phase shifts in 
the two spin channels, rjs + i/ 2 {u) and rj S _ 1/2 (w), from the 
energy levels of the low-lying single particle excitations in 
the approach to the fixed point (N — > oo). The relation 
of the phase shift to the NRG flows is 

E p , 1 (N)=E* 1 (l-±r 1 (N)) (17) 

with El denoting the one-particle excitation energy at 
the fixed point. Translating the N dependence in a lu de- 
pendence and making use of the fact that E^l (N) — > E{ 
very rapidly, we can use Eq. J7J) and Eq. (JSJ) to calculate 
the phase shift. For the 5 + 1/2 channel we obtain 

to+i/aM = -^ttSJH (18) 

and similarly for the 5 — 1/2 

Vs-x/2^) = ^(S+l)J(u J ). (19) 

The numerical value of the proportionality factor is 
®0i/E* = 0-49 for a A = 1.8. In appendix lAl we con- 
sider explicitly the exchange scattering of a single quasi- 
particle by an interaction </, which should asymptotically 
describe the behaviour near the fixed point. The result- 
ing phase shifts are given in Eq. (| A 18(1 and can be seen 
to give the same as those estimated from the level shifts 
for w -> (J(uj) -> 0). 

B. Antiferromagnetic coupling (J > 0) 

For antiferromagnetic coupling, a partial screening of 
the spin 5 takes place yielding a ground state with a re- 
duced spin 5 — 1/2. This screening introduces a phase 



shift of 7r/2, and therefore, in order to facilitate a com- 
parison with the ferromagnetic case, we consider only odd 
iterations here. 

We assume that the low-lying excitations from the 
ground state can again be expressed using Eq. Q and 
calculate J(N) using the flows of the NRG. Figure |3| 
shows the flow of 

E Ptl (N, 5) - E pA (N, S-l) = J(iV)(25) a 2 1 A( JV - 1 )/ 2 

multiplied by N. This corresponds to the quantity ana- 
lyzed in the ferromagnetic case as the antiferromagnetic 
model has a ground state with S replaced with 5" — 1/2. 

Looking first at the curves for smaller S, we see two 
distinct regions. Initially, J(N) x N decreases linearly 
with N. For larger values of N we see the flows increasing 
relatively rapidly. For higher values of S, the flow is 
monotonically decreasing as in the ferromagnetic case. 
This suggests that the low-energy fixed point is the same 
as in the ferromagnetic case but with an additional phase 
shift of 7r/2 and a spin S — 1/2 reduced by 1/2. Also 
the approach to the fixed point is formally identical, see 
Eq. with S S- 1/2. 

As can be seen in Fig. [3] and the approach to the 
fixed point can vary depending on J and S. For Tq/D > 1 
(Tq/D < 1), the fixed point is approached from be- 
low (above). Looking at Fig. we see that for a fixed 
5 = 3/2 there is a crossover from a non-monotonic flow 
for smaller values of J to a monotonically decreasing flow 
for higher values. In the latter case, the flows are com- 
pletely similar to those with a ferromagnetic bare J < 0. 
Therefore it makes sense to use Eq. (|14|) to define a ef- 
fective bare J e g for the antiferromagnetic case. We will 
see below that J e g is usually ferromagnetic. Then, using 
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FIG. 5: The effective ferromagnetic coupling J a g as a function 
of the bare antiferromagnetic J for S = 1 (x), S = 3/2 (o), 
and S = 2(D). 



Eq. (fTTTf) . the coupling reads 

3{N) '- N-B(J )MJ) - < 21 » 

The parameter J c g can be interpreted as an effective 
coupling between the residual spin and the conduction 
electrons, provided we are in the strong coupling regime. 
From Table [I] we infer that in this regime a reasonable 
approximation to chose a constant B. Therefore we use 
a fixed value B — 1.62 in order to fit Eq. to the NRG 
flows (Eq. CT). 

The fits yield the same value of A ss 1.57 as in the 
ferromagnetic case for all S and J. The results for J c s 
as a function of J are plotted in Fig. [SJ We observe that 
JcS — ► as J — > oo and that J e ff depends strongly on 
5. For smaller values of J than plotted in Fig. J e g 
is positive and hence the coupling is antiferromagnetic. 
It should be pointed out, however, that J(N) is always 
ferromagnetic. This can also be seen in Fig. 0] where all 
energy differences are negative. 

As in the ferromagnetic case, we can derive the phase 
shifts in the two spin channels. These evaluate to 

Vs{") = \-^{S-\) J» (22) 

in the high-spin case and similarly 

ris-i{u) = \ + ^{S+\)J{u) (23) 



III. GREEN'S FUNCTIONS AND T MATRICES 

For the calculation of dynamic quantities, we apply 
the NRG with a discretisation parameter A = 1.8 and 
retain up to 1600 states. We focus, first of all, on the 
one-electron Green's functions and t-matrices. To calcu- 
late the t-matrix for the s-d model, we take the equa- 
tion of motion for the single-electron Green's function 
Gk-\,wi(u>) = ((c fc | : c\, |)) and obtain the relation 

Gfe Tife ' T (cj) = 1 G t (uj) , (24) 

where G t (to) is given by 

G t (u) = ((S-d L + : S+d\ + S z d\)) . (25) 

Hence in this case the on-shell t-matrb*2& can be ex- 
pressed in terms of the Green's function Gt(u>) via 

tk,k'(uj) = J 2 a k a k > G t (uj) ■ (26) 

We can find a relation between Gt(u>) and the local elec- 
tron Green's function Gd{oj) = ((d a : d' a J) on the d site. 
For this we multiply Eq. (|24|l by a k oik' and sum over all 
k and k'. We find 

G„M = G^V) + J 2 G< 0) M G t (w) Gi 0) (c) , (27) 

with the non-interacting (i-site Green's function 

GJV)=£— ■ (28) 
, w — e k 

k 

In the wide band limit we can take G? (ui) = —iirpo for 
u> <C D, and in this case we find 

G d (cj) = -iirpo - Tr 2 p 2 J 2 G t (Lu) . (29) 

From this we deduce for the spectral densities 

Pd{u) = Po[l - TT 2 paJ 2 p t (uj)] , (30) 

where pt{u) = — Im Gt(u>)/ir. 

We know from the results for the Anderson model (see 
appendix [BJ) that pd(0) = 0. This must apply to the s- 
d model with antiferromagnetic coupling with spin S = 
1/2, because it is equivalent to the Kondo limit of the 
Anderson model. By inserting pd(0) — into Eq. I|30|) 
we find 

Pt(°) = ~T~n ■ ( 31 ) 



in the low-spin case. The former phase shift agrees 
with that calculated from the spinon density of states 
inRef.lU 



We conjecture that this result applies quite generally for 
antiferromagnetic coupling and all values of S. This is 
confirmed in our numerical results given below. 
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FIG. 6: Spectral densities p/ of the Anderson model and p t of 
the Kondo model (top) as well as pd (bottom) for both model; 
A resonance in pf and pt corresponds to an antiresonance 
in pd- 



A. Anderson model and Kondo model 

In this subsection we look at the relation between the 
Kondo model (in our notation: s-d model with 5 = 1/2 
and antiferromagnetic coupling) and the Anderson im- 
purity model. In the latter model the spin S = 1/2 
is replaced with an /-level that couples to the con- 
duction band at the site d with hybridization V, see 
Fig. ^ On this /-site the electrons interact with a Hub- 
bard repulsion U. The Anderson model and the Kondo 
model are known to be related by the Schrieffer- Wolff 
transformation 29 which yields J = AV 2 /U. 

Another way of relating these two models is to look 
at the singlet-triplet splitting at the isolated impurity in 
both cases. For symmetric parameters of the Anderson 
model this yields 



By expanding the square root to first order in V 2 /U 2 one 
recovers the Schrieffer- Wolff result. 

Figure compares the spectra of the two models. The 
top panel shows pt of the Kondo model along with the 
spectrum pf of the local /-site Green's function for two 
different values of U. Dashed lines refer to a value of V 
as calculated from the Schrieffer- Wolff transformation, 
whereas dot-dashed lines are for a V as calculated from 
Eq. (|32|l . One observes that in all cases the central reso- 
nance of the Kondo model is well approximated by pf of 
the Anderson model. V from the singlet-triplet splitting, 
however, yields a slightly better approximation to the 
Kondo model. Obviously, p t of the Kondo model does 
not show the charge excitations present in the Anderson 
model's /?/. 

The same can be said for the anti-resonance at the 
d-site shown in the bottom panel of Fig. This anti- 
resonance is expected from Eq. H3U|) . Again the singlet- 
triplet splitting yields a somewhat better approxima- 
tion to the anti-resonance of the Kondo model than the 
Schrieffer- Wolff result. 



B. Ferromagnetic Coupling 

The case of a ferromagnetic coupling to a spin S =1/2 
is completely different from the antiferromagnetic one, 
since there is no tendency to screen the impurity spin. 
There is no corresponding Anderson model. 

Qualitatively, the shape of the spectra is reversed, as 
can be seen in Fig. [7] For ferromagnetic coupling, p t has 
an anti-resonance whereas pd a resonance whose height 
is independents^ , of J. This independence stems from 
the fact that the anti-resonance goes straight down to 
Pt(0) = 0, in which case Eq. lp?U)) implies pd(0) = pa- 
The low-energy behaviour of pt(w) can be fitted to the 
function 



which vanishes at u> = 0. This point is illustrated in 
Fig. [S] where the rescaled spectra for various coupling 
strengths are plotted on a logarithmic low-energy scale. 
Equation (|33|) also reflects the singular approach of Pt(u) 
to lo = as the derivative diverges at this point. This 
singular behaviour is also seen in the cusp-shaped peak 
of pd at lo = (see inset to Fig. |TJ . 

As is shown in Fig. [51 To seems to be independent of 
S and can be fitted with the formula 

T (J) = D V^W"e _1/(2poJ) (34) 

with D = 2.55 and po = 0.47. This energy scale agrees 
remarkably well with the energy scale To obtained from 
the fixed point analysis, see Eq. In the ferromag- 

netic case, T is a large energy scale and greater than the 
bandwidth D. 
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FIG. 7: Spectra p t (top) and pa (bottom) for S = 1/2 and 
various values of J for ferromagnetic coupling. The anti- 
resonance in pd corresponds to the resonance in pt. 



C. Antiferromagnetic Coupling 

We turn to the spectra of the underscreened s-d mod- 
els with antiferromagnetic coupling to a spin S > 1/2. 
Here, a resonance at u> = is found in Gt(uj) and an 
antiresonance in Gd(u)). These are plotted in Fig. 1101 for 
5 = 1 and various values of J. The features at \u>\ w 1 
are due to the band edge and pt(0) « l/(J 2 7r 2 po) as one 
would expect from Eq. (|31[1 . The gap in Pd{u) at lo = 
remains for all J < 0, and pd(0) = however small J in 
contrast to the situation J = 0, where there is no gap 
and p d (uj) = po- 

Figure 1111 shows pt(w) for various values of S and 
a fixed J. The band edge features become more pro- 
nounced with increasing S, and pt(0) is independent of S. 
A magnified picture of the low-energy behaviour is shown 
in the inset. On this scale, one clearly sees the qualita- 
tive difference between the fully screened S — 1/2 case 
(quadratic behaviour) and the underscreened S > 1/2 
case (cusps). In the latter case, a good fit to the low- 
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FIG. 8: Low-energy behaviour of the rescaled spectral func- 
tion p t as a function of the reduced energy u>/Tq for FM cou- 
pling to a spin S = 1/2. After rescaling, the spectra corre- 
sponding to all coupling strengths fall on one curve that is 
well fitted by p t (uj) = 6/ ln 2 (T /oj). The range of fitting is 
always the interval [10 _8 ,0.5] independent of J which takes 
the values -0.1 (x), -0.15 (o) to -0.45 (<). 



energy behaviour of the spectra pt(ui) can be obtained 
with 



p t {ui) ~ a 



ln>/T ) 



(35) 



The parameter a is independent of S and approximately 
given by a = 1/ (J 2 ir 2 po) as inferred from Eq. I|30|l for a 
vanishing pd(0). In Fig. El we plot the rescaled spectra 
on a logarithmic scale. On low energy scales, all spectra 
collapse onto a single curve showing (i) the validity of 
Eq. H35JI . and (ii) that b is an implicit function of To- The 
resulting energy scale To is shown in Fig. 1 131 as a function 
of J for various values of S. The dashed lines indicate 
To as obtained from the flow diagrams. We observe that 
the agreement of the two approaches to extract To is not 
as good as in the ferromagnetic case. 

The reason for this disagreement is that these two Td's 
are obtained from the behaviour on low but rather differ- 
ent energy scales. There is no a priori reason for them to 
coincided The antiferromagnetic case is more complex 
than the ferromagnetic. A higher energy scale (Kondo 
temperature Tk, independent of S) is associated with 
the partial screening of the original spin. By contrast, on 
the lowest energies the residual ferromagnetic coupling 
leads to another energy scale Tq depending on S, as dis- 
cussed in Sec. [H] A reasonable explanation of the fact 
that T obtained from the spectra differs from the one 
from the fixed-point analysis is that the former is still 
affected by the crossover between these regions. 



8 




0.1 0.2 0.3 0.4 0.5 



Ml 

FIG. 9: Energy scale To as a function of the ferromagnetic 
coupling strength J for different S. Estimates of To from the 
fit of p t (u) = b/\n 2 (io/T Q ) over the interval [1CT 8 , 10" 2 ] are 
shown as (o) for S = 1/2 and as (x) for S = 1. Estimates 
of To from the flows through Eq. 11311 are denoted by o for 
5 = 1/2 and * for S = 1. The full line joining the points is 
given by T (J) = D y/2p^J exp(-l/(2po J)) with D = 2.55 
and po = 0.47. In the inset the parameter b is plotted as a 
function of J. Here the data for S = 1 are roughly 3 x those 
for 5" = 1/2. 



IV. SCATTERING CROSS SECTIONS 

The on-shell t-matrix tk,k' directly yields the elastic 
cross section and can be related to the total scattering 
cross section via the optical theorem. Following Zarand 
et a/—, we can write the scattering cross sections as 




CO 



FIG. 10: Spectral densities p t (top) and pd (bottom) of the 
underscreened s-d model with S = 1 resonance and antireso- 
nance; former normalized by J 2 iv 2 po 



ftot(fc) = Imife,fc 

' d 3 k' (36) 
2ir6(e k > - e k ) \t k ,. k \ 2 , 



tT c i(fc) = 



1 



v F 



(2tt) 3 



where vf denotes the Fermi velocity. The inelastic cross 
section is then given by the difference ai ne \(k) = <Jtot(k) — 
a e i(k). By evaluating these equations for a flat density of 
states po and ty,k = J 2 ct k ce k > Gt(e k ) with a k depending 
on k only through e k , we obtain 



ototM 



2J 2 \a k \ 

VF 



■lmG t (u) 



2WP ° N J*\a k \*\G t (u,)\* 



(37) 




This simplifies further if we assume purely isotropic scat- 
tering a 2 , = 1/N, and we find 

otot M = -voPon J 2 Im G t (uj) (38) 
CTciM =a pW J 4 \G t (iu)\ 2 (39) 



FIG. 11: Spectral densities pt of the antiferromagnetically 
coupled s-d model with S = 1/2 (dotted), S = 1 (dashed), 
5 = 3/2, S = 2 (dot-dashed), and S = 5/2. The inset shows 
the resonance on a lower energy scale, where the spectrum of 
the fully screened model (dotted) is almost flat. The value at 
lj — is approximately the same for all values of S. 



9 



0.95 



0.9 



„- 0.85 



0.8 



0.75 



0.7 L 



10 



10 



CO / co„ 



10 



10 



FIG. 12: Rescaled spectral densities pt of the AFM s-d model 
with S — 1 (top) for antiferromagnetic couplings J — 0.1(x), 
J = 0.15(o), . . ., J — 0.45. One sees that after rescaling 
the spectrum p and the frequency uj, all points fall on the 
same line. For better visibility, the inset shows the same data 
slightly shifted for each J. Full lines are least square fits to 
a — bj ln 2 (Tb/u;). The fits have been done on the interval 
[10 -8 ,u;o], where loq is taken such that p t (oJo)/ Pt(0) « 0.9. 




FIG. 13: Energy scale To as a function of the antiferro- 
magnetic coupling J for S = 1 (diamonds) and S — 3/2 
(squares), and S = 2 (pentagrams). The full lines show fits 
to T (J) = D \bJ\ n exp(6/J) with different values of D, b and 
n. The dashed lines with the somewhat smaller symbols indi- 
cate To as extracted from the flow diagrams. The dot-dashed 
line shows Tk = D v / 2poJexp(— l/(2po J)) with D = 1 and 
po = 0.5. 




FIG. 14: Elastic (dashed), inelastic (dot-dashed) and total 
(full) cross section for the fully screened model (left) and un- 
derscreened model (S = 1, right) for J = 0.2. 




FIG. 15: Elastic (dashed), inelastic (dot-dashed) and total 
(full) cross section for the ferromagnetic model J = —0.2 The 
left-hand (right-hand) plot is for S = 1/2 (S = 1). The dot- 
dashed and the full lines virtually lie on top of each other 
indicating that the scattering is predominantly inelastic. 



Therefore, on the Fermi level 



<7tot(0) = <7 c i(0) 



ttpqvfN 



00 



(40) 



This shows that the inelastic scattering cross section van- 
ishes at the Fermi level, as it should for a Fermi liquid. 
For the regular Fermi liquid (S = 1/2, Kondo model), 
the inelastic scattering vanishes quadratically, and ther- 
fore on the lowest energy scale tTtot(w) = a e \(u>). In con- 
trast to that, for the underscreened model (S = 1), the 
inelastic scattering falls off much more slowly. 

In the ferromagnetic case, we find that Gt{oj) — > as 
Cij — > 0, which implies that the total cross section and 
thus all scattering vanishes at the Fermi level. Moreover, 
since cr e i(w) involves the square of Gt(tu) <C 1 for uj — * 
we expect (t c \(uj) <C (Ttot(^) ~ Gt{oj) which is a formal 
explanation of the predominance of inelastic scattering 
in this case. 



with (Tq 1 — irpoVpN/2. These cross sections are shown in 
Fig. 1141 for an antiferromagnetic coupling and in Fig. 1151 
for the ferromagnetic case. 

In the antiferromagnetic case, the asymptotic be- 
haviour of G t {oj) as uj — > is G t {uj) — > —i/(poTrJ 2 ). 



V. SPIN SUSCEPTIBILITIES 

From the NRG we also calculate the impurity spin sus- 
ceptilibities x( w ) both for ferromagnetic and antiferro- 
magnetic coupling. The low energy behaviour of x( w ) 
can be obtained from the equations of motion. As is 
shown in appendix [UJ the low-energy behaviour is also 
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FIG. 16: Low-energy spectra of the impurity spin susceptibil- 
ity multiplied by ui for ferromagnetic coupling to S = 1/2 
(dark symbols) and S — 1 (light symbols). The energies 
are rescaled by T x which is fitted from — ui Imx(u)/ir = 
a I \v?(u)/T x ) with a logarithmic least squares fit. Eg. 1431 The 
ranges of the fits are [10~ 8 : 2 x 10" 2 ] for 5=1/2 and [1CT 8 : 
2x 1(T 3 ] for S = 1. J ranges from J = -0.05(*), -0.1(x), . . . 
to J = -0.45(>). 



FIG. 18: Low-energy spectra of the impurity spin susceptibil- 
ity multiplied by uj for antiferromagnetic underscreened mod- 
els S — 1 (dark symbols) and 5 = 3/2 (light red symbols) and 
5 = 2 (light grey symbols). The energies are rescaled by T x 
which is fitted from -wlmx(w)/ir = a/ \v?(lo /T x ) with a log- 
arithmic least squares fit Eq. il l-'il . The ranges of the fits vary 
significantly as a function of J to fit the low-energy behaviour. 
J ranges from J = 0.05(*), 0.1(x), ... to J = 0.45(>). 
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FIG. 17: Energy scale T x as a function of J for ferromagnetic 
coupling of a spin 5 = l/2(o) and 5 = l(x). The broken line 
is a fit to Eq. <E2J| with a = 4.7 and b = 0.46. The full line 
is a better fit to Eq. 1421 with (bJ)~ 1//5 instead of the square 
root and D = 2.55 and b = 0.50. 



singular. It goes asymptotically for u> 
1. 5(5 + 1) 



as 



Im((5 z : S z )) 



3wln (u)/T x ) 



(41) 



1 



ln(ln(a;/r x )) 
H">/T x ) 



o 



(ln(ln(o;/r x )) ; 
ln 2 (o;/T x ) 



where we have introduced an energy scale T x which we 
expect to vary as 



T x = Dy/2po\J\, 



-l/(2Jp ) 



(42) 



This leads to the following ansatz for the low-energy be- 
haviour of the spectrum of x( aJ ) : 



Imx(w) = 5 

win (u>/T x ) 



(43) 



As can be seen from Figs. 1161 and 1181 this formula de- 
scribes the low-energy susceptibilities in both the ferro- 
magnetic and the antiferromagnetic cases very well. In 
the ferromagnetic case, a w 5(5 + l)/3 very well with 
an error of less than 1%. In the antiferromagnetic case, 
5 needs to be replaced by 5 — 1/2. However, here the 
agreement is less accurate with an error of the order of 
10%. 

The characteristic energy scale T x is plotted in Figs. 1171 
and EH We find that T x ps 2 x T both for ferromagnetic 
coupling and also for antiferromagnetic coupling, with To 
from the fitting of the Green's functions over the same 
energy scale. 



VI. CONCLUSIONS 

Our NRG calculations give a comprehensive picture of 
the dynamics of the underscreened and the ferromagnet- 
ically coupled impurity s-d models. In the ferromagnetic 
case, the ground state has (25 + l)-fold degeneracy as 
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FIG. 19: Energy scale T x for underscreened s-d models with 
S = 1 (o), S = 3/2 (□), and 5" = 2 (*). The values of T x are 
obtained from low-energy fits to Eq. 1431 . Full lines are fits 
to Eq. 142t with the square root replaced by the exponent 
n = 2.5 for S = 1, n = 5.0 for S = 3/2, and n = 7.5 for 
5 = 2. 



the impurity spin becomes completely decoupled. The 
NRG excitations, however, only slowly approach their 
fixed point values as 1/(N + C(J)) as N — > oo. This be- 
haviour is consistent with a frequency dependent renor- 
malized exchange coupling J{u>) = l/ln(w/To), which is 
a marginally irrelevant operator at this fixed point. The 
energy scale To, which we calculated in terms of the bare 
parameters was found to be independent of the spin S. 
Its value, as well as the T x deduced from the dynamic 
susceptibility, corresponds to the usual expression for Tr-, 
but as J < 0, it is such that Tk ^> D. 

Essentially the same kind of behaviour was found in 
the case of the underscreened antiferromagnetic model. 
After the initial screening of a component of the spin by 
the conduction electrons in the single channel, the ef- 
fective spin value is reduced to S — 1/2. The effective 
exchange coupling J(u>) on the low energy scale was also 
found to be of the form, J(ui) = 1 / \n(uj / T ) , and hence 
is ferromagnetic for u> < Tq. The values of To obtained 
in this case, however, were found to depend on the spin 
value S, and to be consistently smaller than in the fer- 
romagnetic case, such that To ~ D. The value of Tq 
does not coincide with the formula for the antiferromag- 
netic Tk which is independent of S and always smaller 
than our values for T . The latter, however, does agree 
well with the energy scale T x deduced from the dynamic 
susceptibility. 

For large values of the bare antiferromagnetic interac- 
tion J we were able to estimate the value an effective 
bare ferromagnetic interaction J c g which would lead to 
the calculated values of T . From the NRG energy lev- 
els in the approach to the fixed point, the asymptotic 
form of the phase shifts as u> — > due to the elastic scat- 



tering of a single particle excitation in the two channels 
S r ±l/2 by the residual spin S r could be estimated, where 
S r = S in the ferromagnetic case and S r = 5— 1/2 for the 
antiferromagnetic. These were found to be in complete 
agreement with the values from an explicit calculation of 
the scattering using the effective low energy Hamiltonian 
(see Appendix lA"|l . 

The singular form of the effective coupling J(u>) is 
manifest as sharp cusps at the Fermi level in the spec- 
tral densities calculated for the Green's functions Gd(w) 
and Gt(u>). For antiferromagnetic coupling the value of 
Pt(u>) at to — was found not to diverge and to approach 
the value 1/tt 2 pqJ 2 independent of the spin value S, as 
conjectured in Eq. (|31f) . For ferromagnetic coupling the 
value of pt(0) was found to be zero for finite J but equal 
to S(S + l)po for J = 0. This clearly demonstrates that 
J = is a singular point, and the behaviour is discon- 
tinuous as a function of J when J — ► when approached 
from the ferromagnetic as well as from the antiferromag- 
netic side. The approach of Prf(w) and pt(<-o) to their 
values at lu = could be described well by terms of the 
form l/ln 2 (w/T ) in all cases where there is incomplete 
screening of the local spin. The energy scale deduced 
from this fitting broadly agrees with the one from the 
fixed-point analysis. In the case of ferromagnetic cou- 
pling the agreement is rather precise. 

The results of the calculations for the inelastic and elas- 
tic scattering cross-sections as a function of frequency ui 
are consistent with the classification of the low energy be- 
haviour of these systems with unscreened residual spins 
as singular Fermi liquids, as proposed by Mehta et aiSS. 
The inelastic scattering goes to zero as co — * as to be ex- 
pected for a Fermi liquid, in contrast to non-Fermi liquid 
behaviour which is characterized by finite inelastic scat- 
tering at u> = Q. The approach of the inelastic scattering 
component to zero in the underscreened cases, however, 
was found to be anomalously slow compared to the corre- 
sponding results for the fully screened antiferromagnetic 
model with S = 1/2. The singular nature of the low en- 
ergy scattering is also evident in the contrasting results 
for the spectral densities of pt{u>) as u> — > in the antifer- 
romagnetic case for S = 1/2 and S > 1/2, shown in the 
inset of Fig. ^2 I n the case of ferromagnetic coupling, in- 
elastic scattering gives the dominant contribution to the 
total scattering cross section. 

To summarize: our results are in broad agreement with 
the conclusions of previous studies of the dynamics of 
underscreened s-d modelsi2*i2*i2i2£. However, as well as 
fully analysing the fixed point behaviour we have been 
able to make precise predictions for the dynamics of a 
range of physical response functions, and have also been 
able to calculate the explicitly the renormalized energy 
scales T), T x etc in terms of the parameters of the bare 
model. 
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33 We could calculate the spectra down to cj = 10 -8 but 
the flows to 10 -60 . Therefore, the fitting of the Green's 
functions is done over an energy range of [10 -8 , 10 -5 ] and 
that for the flows on [10 -60 , 10" 40 ]. 



APPENDIX A: EXCHANGE SCATTERING IN 
THE ONE-ELECTRON CASE 

Consider a system of free electrons that couple to 
a spin 5 with a local exchange scattering term. The 
Hamiltonian is given by the sum of Eq. (JJJ and Eq. J5J • 
We calculate the one-electron states using the basis set 
{\k,a : S,M S ) = c\ \0) \S, M s )}. The Hamiltonian acts 
on these states as 

H\k, |: 5, M s ) = e k \k, |: 5, M s ) + a k jj^ a k' 

k' 

x (y{S - M S )(S + M s + 1) \k',i: 5, M s + 1) ( A1 ) 

+ Ms\k',V-S,M a )] 
and 

H\k, j: 5, M s + 1) = e k \k, j: 5, M s + 1) + a k Ofc' 

k' 

x (y(S-M s )(S + M s + l) \k', T: 5, M s ) 

-(Ms + l)\k',i:S,M s + l)) . (A2) 

Solving the eigenvalue equation we find the Ms depen- 
dence cancels out, and the equation factorizes so that we 
get two independent equations for the new energy lev- 
els. With the free ci-site green's function as defined in 
Eq. we obtain 



1 = JSG { °\u) 



(A3) 



which corresponds to scattered states in the S+l/2 chan- 
nel and 



1 = -J(5 +!)£<%) 



(A4) 



which correspond to scattering in the S — 1/2 channel. 

The one electron eigenstates in S+l/2 channel are 
built up from the states, 



|fc,M s + l/2) = 



5 + M s + 1 



25+1 



1/2 



\k,V.S,M s ) 



S-M s 
25 + 1 



1/2 



(A5) 



\k,l:S,M s + l) 
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and those for the 5—1/2 channel from the states, 



\k, M s + 1/2) = - 



5 -M s 
25+1 



1/2 



\k,t-S,M s ) 



S + Ms + l 
25+1 



1/2 



\k,i: S,M S + 1) 



(A6) 

We can define corresponding creation operators, 
via 



fe,M s + l/2 



4,m s+ i /2 = v 25 A +r 1 /(5 ' Ms)(5+)(Ms+s)c M 



S-M s 
25+1 



(A7) 

/(5,Af s + l)(5+)^+ s+1 ) Cfc4 



where Mg can only take values 5, 5 — 1,...,— 5 — 1 as 
(5+) 2S + 1 =0,and B[ Ms+l/2 via 



D t,M s + l/2 



Vf5TT /(5 ' Ms)(5+)(Ms+S)c ^ 



(A8) 



+ f\slV f ^ Ms + 1)(* + ) (Ms+5+1) cm 

where M$ can take values 5 — 1,...,— 5. These 
operators operate on the vacuum state |0)|5, — 5), 
and /(5, Ms) is a normalization factor such that 
/(5,M s )(5+)( Ms+s )|5,-5) = \S,M S ). 

Confining our attention to the system with no more 
than one electron we find equations for the Green's func- 

ti0nS G fc,fc',M s +l/2M = {{ A k,M s +l/2 : A l',M s +l/2» aild 
^k,k' jMs+y^) = (( B k,M s + l/2 : B k' ,Ms + l/2)) aS 

(w - e k ) G^ k , Ms+1/2 (uj) = 5 klk ' (A9) 

+ a k JS^2 a k" Gk",k',Ms + l/2 M 



and 



(uj - e k ) G^ k , Ms+1/2 {u) = S k ,k> (A10) 
-a k J(S + l)^ afc " G fe",fc',M s +i/ 2 ( w ) ■ 



Solving these equations we obtain 

Sk,k' 



Gk,k',M s + l/2( U} ) 



and 



Cfc,fc',M s + l/ 2 ( a; ) 



w - ey 
Qfc J5 a k > 

8k, k< 



(All) 



uj - ek' 
a k J(5 + l) 



(A12) 

OL k i 



ui — 



£ fel + J{S+l)G [ °\uj)u-tk> 



This enables us to calculate also the Green's function 
G t (u>) as defined in Eq. 1251) for the one-electron case. In 
terms of the creation and annihilation operators defined 
in equations (|A7I) and (|A8|) . we find 



1 T Y V25TT 



(A13) 



1/2 



; {5V5 + M s + l4 Ms+1/ 
-(5 + l)VS^Ms4 iMs+1/2 }. 



The Green's function Gt(uj) for this one electron situation 
is then given by 



OtkUk' 



^ 2(25+l) 2 

k,k',M s y ' 



(A14) 



{ 



5 2 (5 + M s + l)G^„ !Ms+1/2 H 
-(5 + l) 2 (5-A/ s )G^, !Ms+1/2 H} 



We find 



E a 1 \ G^ ^ (uj) 
aka k iG kk , Ms+l/2 (uj) = - Tan(f))i 

k,k' 



l-JSG$>(u) 



(A15) 



and 



a kak'G^ k , Ms+1/2 {uj) = 



k ,k> 1 + J(S+1)G { °\uj) 

(A16) 

We substitute these into Eq. I|A14(1 and perform the 
sums over Ms- This sum has to be done carefully be- 
cause the for the first term the sum is over the values 
— 5 — 1, — 5, . . . , 5, whereas for the second term it is over 
-5, -5 + 1, . . . , 5 - 1. The result is 



5(5 + l)Gi 0) M 



[l-JSGf\uj) 1 + J(5 + 1)G^HJ 

In our many-electron case this theory is only directly ap- 
plicable in the immediate vicinity of the fixed point with 
one quasiparticle excited from the interacting ground 
state and J interpreted as the renormalized coupling J. 
The phase shifts for the quasiparticle scattering are then 
given by 



(A17) 



(0), 



7 ?S+i/ 2 ( w ) = -arctan(7rpoM</(w)5) , 
175-1/2 (w) = arctan (ttp (uj)J(u))(S + 1)) 



(A18) 



with J (uj) = 1/ ln(cj/Tb). In the limit id — > 0, the phase 
shifts we calculate in this way agree with the estimates 
Eq. O and Eq. O obtained from the NRG flows. 



14 



APPENDIX B: GREEN'S FUNCTIONS AND 
T-MATRICES FOR THE ANDERSON MODEL 

Let the impurity site be labelled by / and next site 
it is coupled to on the chain be labelled d, as indicated 
in Fig. EJ Let d a = E k V k c k /V where \V\ 2 = £ fc \V k \ 2 . 
The equation of motion for the Green's function G kk > (lo) 
yields the relation 



G kk > (u>) 



5 kk > 



cu - e k uj - e k 



■ G /M - — ~J 



The t-matrix t kk i (w) is defined by the equation, 

Skk' 



G kk > (ui) 



lo - e k u> - e k 



1 , 1 

■ t kk >(U)) 



(Bl) 



(B2) 



and hence for the Anderson model we have 

t kkl {uj) = V 2 a k G s {uo)a k , . (B3) 

We can derive an expression for Gd(<jj) from G kk '(u) 
by multiplying by V k , V k > , and summing over k and k! . 
This gives 

G d {u) = G< 0) ( W ) + V 2 G$\ U ) G f {u>) V 2 G<%) . (B4) 
In the wide band limit V 2 G^ (0) = — iA and hence, 



1 



G <*M = ^ [-iA- A 2 G f {u>)] 



(B5) 



We can deduce the spectral density of Gd(u>) from 
1. 



Pd{u) 



-ImG / H=p [l-^Ap/H], (B6) 



where we have used the result A = 7rpo|V| 2 ! where po 
is the free conduction density of states at the Fermi 
level. It follows from this expression that for the sym- 
metric model, where p/(0) = 1/ttA, the spectral density 
of Gd vanishes, pd(0) = 0. Therefore, we find an anti- 
resonance. 



APPENDIX C: DYNAMIC SPIN 
SUSCEPTIBILITY 

We derive equations of motion for the local spin-spin 
double-time Green's function {(S + (t) : S~(t'))) to second 
order in the coupling J. For this we need the evaluation 
of the higher order Green's function to zero order only. 
We find from the equations of motion 



oj 2 ((S+ :S-)) U = 



8J 2 S(S + 1) 



((si ■■ Sj)) u , (CI) 



where we have used the fact the ((S + : S~)) = 2((S Z : 
S z )), due to the rotational symmetry in zero magnetic 
field. The spin operator Sj" for the first site on the chain 
in terms of the conduction electron states is given by 
= d^-d^ = J2 k k' a fc Q; fe' c fcT c fe'i' w ith a similar expres- 
sion for S7. Evaluation of ({St : Sj)) to zero order in 
J gives 



«S* ■ S d)) ^ /J Po(e) PO {e') ( 



n{f(e)-m 



(j — e + e 



de de' 



( C2 ) 

where po(u)) = J2 k a fe^( w — e ) i s the local spectral den- 
sity for the first site on the chain for the non-interacting 
system. From that we deduce that the spectral density 
of ((S + : S~)) to order J 2 is given by 



1 + 8J 2 S(S+l) 
—lm((S + : S )) = —2 

x / po(e)po(e + (j)(f(e)-f(e + u>))de. 



(C3) 



For the limit u> — > we find 
-±^^ {{ S+:S-)) = 8 ^ 0)) y + 1) (C4) 

7T OJ^O 3 

Using the poor man's scaling equation, Eq. (3.51) in 
Ref. l32l with D — au>, the coupling J is renormalized to 
J(u>) given by 



2J Po 



1 



HHlo/T x ) 



o 



1 



ln(w/T x ) 2\n 2 (uj/T x ) \\ii 2 (uj/T x ) 
with the energy scale T x given by 

T x = D\2Jp \ 1/2 e- 1/2Jpa . 



(C5) 
(C6) 



Hence, asymptotically as u> — > 0, 



1) 



--lm((S+ : S-)) = 

7T 3wln 2 (w/T x ) 



1 



ln(ln(c/T x )) 
ln(w/T x ) 



O 



(ln(ln(^/r x )) ; 
ln 2 (u;/T x ) 



(C7) 

For the fully screened Kondo model, this asymptotic form 
should be appropriate only for the regime D > oj > T x — 
Tk- However, for a ferromagnetic coupling and for the 
underscreened models, this asymptotic form should also 
apply as u> — > 0, as confirmed in Sec. |VJ As the asymp- 
totic behaviour of the underscreened models is associated 
with a ferromagnetic fixed point, the low-temperature 
T x is expected to differ from Tk (high temperature) , see 
Sec. El 



